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ABSTRACT 

We apply the matrix model of Kapustin, Willett and Yaakov to compute the 
free energy of Af = 3 Chern-Simons matter theories with D n quivers in the large 
N limit. We conjecture a general expression for the free energy that is explicitly 
invariant under Seiberg duality and show that it can be interpreted as a sum over 
certain graphs known as signed graphs. Through the AdS/CFT correspondence, 
this leads to a prediction for the volume of certain tri-Sasaki Einstein manifolds. 
We also study the unfolding procedure, which relates these D n quivers to A 2n -b 
quivers. Furthermore, we consider the addition of massive fundamental flavor 
fields, verifying that integrating these out decreases the free energy in accordance 
with the F-theorem. 
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I 



1 Introduction 



Supersymmetric localization [1, 2] is a powerful method that makes exact computations in 
superconformal field theories (SCFTs) possible. This procedure reduces the infinite dimen- 
sional path integral to a finite dimensional integral, typically over the Coulomb branch. 
It has recently been used to obtain interesting results for field theories in various number 
of dimensions [2, 3, 4, 5, 6, 7]. In particular, Kapustin, Willett and Yaakov [3] applied 
localization in three dimensions to calculate the exact partition function on S 3 for theories 
with M > 2 supersymmetry. One of the outcomes of these calculations has been the proposal 
[8] that the free energy, defined by 

F = -log|Z 5 3|, (1.1) 

decreases along renormalization group (RG) flows, providing a good measure of the degrees 
of freedom in the field theory. On the other hand, many three-dimensional SCFTs can be 
realized as effective theories of coincident N M2-branes. Thus, localization is a tool that 
can be used to test predictions of the AdS/CFT correspondence [9, 10, 11]. One of the first 
and most remarkable results [12] was the evaluation of the free energy for U(N)k x U(N)^k 
Chern-Simons (CS) theory [13], matching the famous N 3 / 2 scaling of gravitational free energy 
predicted in [14]. 

A larger class of quiver Chern-Simons theories with U (iV)&i x U O^Ofo x ... x [/ (iV)k n gauge 
groups, coupled to bifundamental matter forming a necklace-type quiver were considered in 
[8, 15, 16, 17]. It is believed that the M-theory description of these theories [18] arises as the 
near-horizon limit of a stack of iV M2-branes placed at the tip of a Calabi-Yau cone with 
a tri-Sasaki Einstein base Y. In the large N limit, the gravitational free energy is given by 
[12, 15] 

where Vol(y) is the volume of the compact manifold Y whose geometry depends on the 
quiver data, in particular the CS levels. By evaluating the free energy (1.1) for the necklace 
quivers and matching it with the gravitational energy, an expression for Vol(Y) as a function 
of the CS levels was conjectured in [15]. This was corroborated in [19] by comparison with 
the explicit calculations of the volumes of toric Sasaki-Einstein manifolds [20] (see also [21] 
for a calculation in type-IIB supergravity) . 

These necklace quivers are actually an example of a more general class of quiver theories 
which have a nice large N limit, i.e., long-range forces between eigenvalues in the matrix 
model cancel [22]. In fact, quiver theories for which this happens are in one-to-one corres- 
pondence with the extended ADE Dynkin diagrams with necklace quivers corresponding to 
the A-class. 

In this paper we focus on theories with D n quivers. The relevant tri-Sasaki Einstein 
manifold Y is the base of the hyperkahler cone EI 4n ~ 8 / / /[/(l) n_1 x SU(2) n ~ 3 . By assuming 
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a particular ordering of CS levels, we solve the matrix models for various values of n and 
propose an expression for Vol(y) for arbitrary n. This expression is related to the area of a 
certain polygon as in the case of A-quivers. Then, we propose a general expression given by 
a rational function, which is valid for any ordering of the CS levels and is invariant under 
Seiberg duality. The numerator for such an expression was given in [22]. Here we give the 
denominator as well, leading to the conjecture 

Vol(y) T.n + det(ai...a w ) 2 ]\ n a=l \a a ■ p\ 

VolCT "8(n-2)(ELxbn|)nL 1 [E6=i(k-p 6 | + ba+P6|)-4b a |]' ( • } 

where 1Z + is an n-subset of positive roots a a of D n and the CS levels are k a = a a ■ p. We 
show that the numerator of this expression can be expressed as a sum over certain graphs 
known as signed graphs. Using a generalized matrix-tree formula, we show that (1.3) reduces 
to the polygon formula for a particular ordering of the CS levels. Although we do not 
discuss exceptional quivers in detail, we give the free energy for Eq, E$ in Appendix C for 
completeness. 

This paper is organized as follows: In Section 2 we review the matrix model, the large TV 
limit and the ADE classification. In Section 3 we explicitly solve the D n matrix models for 
n = 5 , 6 and conjecture an expression for arbitrary n. In Section 4 we propose the general 
formula and the relation to signed graphs. In Section 5 we add fundamental flavor fields and 
check the F-theorem. In Section 6 we apply the procedure known as unfolding, which relates 
the free energy for .D-quivers to that of A-quivers for a particular choice of CS levels. We 
conclude with a summary and discussion of open problems. 



2 ADE Matrix Models 



We will consider quiver Chern-Simons gauge theories involving products of unitary groups 
only, i.e., G = ® a U(n a N), coupled to bifundamental chiral superfields (A a , B a ). According 
to [3], the partition function of these theories on 5* 3 is localized on configurations where the 
auxiliary scalar fields a a in the M = 2 vector multiplets are constant N x N matrices. Thus, 
evaluating the free energy amounts to solving a matrix model. 



Matrix Model 



We denote the eigenvalues of a a in each vector multiplet by A 0) j, i = 1, N a . The partition 
function is then given by [3] 

Z = j {^\ dX -^j L v ({X a ,i})L m ({Ki}) = j (jJd\ a ,ijexp[-F({\ aii })}, (2.1) 
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where the contribution from vector multiplets is 



Lv = II aTi ( II 2 sin Hn(K,i ~ K,j)) J exp | in £ k a X 2 ay 

a=l a \i>j / \ a,j 



and from matter multiplets is 



L m ~ n n^—Tr-fi — r Ti n (n 



(a,6)eS 



^- 2 cosh[7r(A aj j — A&j)] ^ 2 cosh[7rA C) , 



The first product in L m is due to bifundamental fields while the second one is due to funda- 
mental flavor fields, where n{. is the number of pairs of flavor fields at the node labeled by 
the index c. 



Large N Limit and ADE Classification 

Following [15, 22], we assume that the eigenvalue distribution becomes dense in the large N 
limit, i.e., A tt) , — > X a (x) with a certain density p(x). In this limit the free energy becomes 
a 1-dimensional integral which we evaluate by saddle point approximation. We also assume 
that the eigenvalue distribution for a node with iV a = n a N is given by a collection of n a 
curves in the complex plane labeled by X a ,i(x) with I — 1, n a and write the ansatz 

X a ,i(x) = N a x + iy a j{x). (2.2) 

The density p(x) is assumed to be normalized, i.e., 

dxp(x) = 1 , (2.3) 



which will be imposed through a Lagrange multiplier p,. As explained in [22], the leading order 
in N in the saddle point equation is proportional to the combination 2n a — Ylb\(a fc)e£ n b- The 
requirement that this term vanishes is equivalent to the quiver being in correspondence with 
simply laced extended Dynkin diagrams, leading to the ADE classification. To next order 
in N, the saddle point equation contains a tree-level contribution and a 1-loop contribution. 
Assuming that ^ a n a k a = 0, the requirement that these two contributions are balanced leads 
to a = 1/2, which is ultimately responsible for the N 3 ^ 2 scaling of the free energy 1 . Finally, 
the Lagrangian to be extremized reads 



F = N 3/2 J p(x) 



d n a n a 



"a / \ I a a 

nn F \x\ + 2nxJ2Y, W*) + 4^ £ £ E ar § (e^'^-^f 



An v 

a 1=1 \a=l 1=1 J=l 



n a n b 

2 



J2J2^( e 



^i{Va,I-Vb,j) 

(a,b)&E 1=1 J=l 



(2.4) 



1 Alternatively, one can assume that ^2n a k a ^ 0, and choose a — 1/3, which leads to a massive IIA 
supergravity dual [8]. We will not consider this case here. 
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where rip = Yl a n a n a- Evaluating the free energy on-shell gives 

F = ^—/x, (2.5) 

which can be understood as a virial theorem [19]. Thus, the free energy is determined by /i, 
which in turn is determined as a function of the CS levels from the normalization condition 
(2.3). Note that from (1.2) and (2.5), it follows that 

M = l (26) 

As mentioned earlier, theories with A m _i quiver diagrams have been extensively studied. 
Here, we wish to study theories with D n quivers as the one shown in Fig. 1. For now we will 
set nl = 0, but we will reintroduce flavors in Section 5. 



Figure 1: D n quiver diagram. Each node 'a' corresponds to a U(n a N) gauge group with CS 
level k a , where n a is the node's comark and we assume that Y2 a n a-^a = 0. 

It is convenient to relate the CS level kt a ) at each node to a root a a , by introducing a vector 
p and writing kr a \ = a a ■ p . This way, the condition Y2 a n a,k a — is satisfied automatically. 
Choosing a basis for the roots of D n (see Appendix A for conventions), we have 

h = -(pi+P2), k 2 =pi-p 2 , k 3 =p n _ 1 -p n , k 4 = p n _ 1 +p n , 

ki = Pi_ 3 - Pi-2 ; i = 5,...,n + l. (2.7) 

In the next two sections we will solve the matrix models for various D n quivers and 
conjecture a general volume formula for arbitrary n. 



3 Solving the Matrix Models 

Here we describe the saddle point evaluation of the free energy (2.4). We show in detail 
the solution for n — 5, state the result for n — 6, and propose a general expression that we 
have checked for n = 7, 10. Finally, we will relate this expression to the area of a certain 
polygon. 
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3.1 Explicit Solutions 



Extremizing (2.4) (with respect to y a j and p) requires an assumption on the branch of the 
arg functions. We will always take the principle value and therefore we assume that 

\y a ,i - y a ,j\ < i ; \y a ,i - Vb,j\ < x > if (a, 6) e £7. (3.1) 

Based on numerical results [15, 22], we assume that the n a curves for a given node initially 
coincide, i.e., \y a j — y a ,j\ = 0. Extremizing F under these assumptions, one finds that 
the solution is consistent only in a bounded region away from the origin. This is because 
as \x\ increases, the differences \y a j — y^j\ monotonically increase (or decrease), saturating 
one (or more) of the inequalities assumed in (3.1) at some point. The relation among the 
CS levels determines the sequence in which these inequalities saturate. This saturation will 
be maintained beyond this point, requiring the eigenvalue distribution involved either to 
bifurcate or develop a kink. As an example, consider the first plot in Fig. 2 where we show 
the eigenvalue distributions for the D§ quiver 2 . It consists of seven regions determined by 
saturation of different inequalities. At the end of first region (x = Xi), one can see that 
- 2/5,2 = -1/2 forcing y 5>1 and y 5>2 to bifurcate. 




Figure 2: The eigenvalue distribution y a j(x) (left) for all nodes and density p(x) (right) for 
the D 5 quiver with CS levels: (k 2 , k 3 , fc 4 , k 5 , h§) = (2, 2, 3, 4, 4). 

After a saturation occurs, the total number of independent variables is reduced by one. 
Thus, at this point, we remove one variable from the Lagrangian, revise the inequalities and 
solve the equations of motion again until a new saturation is encountered. This process is 
iterated until all y a 's are related, determining a maximum of {Y^ a n a — 1) regions or until 
the eigenvalue distribution terminates, i.e., p(x) = 0. Once the eigenvalue density p(x) is 
determined in all regions, the value of p (and therefore F) is found from the normalization 
condition (2.3). 

2 We have used the freedom to add an arbitrary function to the y a j to set yi ! i(x) = in the first region 
and we solve explicitly only for x > since the eigenvalue distributions and density are even functions of x. 
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Region 

i 



M 

3(fc2+fc3 + fc4+2fc 5 +2fc 6 ) 

4^ 

6fc^+9fcM+9fc 4 +12fcr,+18A:fi 



_2fj_ 



3(2k 2 +k a +k 4 +2k 5 +2k 6 ) 



2fi 



2k 2 +3(k 3 +k 4 +2k 5 +2k 6 ) 



2ji 



2fc^+3fc3+5fc4+4fc s +6fcfi 



2fi 



2fc 2 +5fc3+3fc4+4fc 5 +6fcB 



2/j 



2fc 2 +3fa +3fc 4 +4fc 5 +6fc fi 



2/1,1 - 2/5,2 = - 
2/5,2 - 2/6,2 = " 



2/2,1 - 2/6,2 = 



2/5,1 - 2/6,2 



2/4,1 - 2/6,2 



2/3,1 - 2/6,2 



2/6,1 - 2/6,2 



1 



\n + \x(ki - k 2 ) 



\x + xki 



H + \x{$k\ — k 3 — 3k 4 — 2fc 6 ) 



2\i + x(2k\ — k 3 — k 4 — k & ) 



Table 1: Key characteristics of the seven regions of the D 5 matrix model: their boundaries, 
the saturated inequalities and the eigenvalue densities, assuming k e > k 5 > k 4 > k 3 > k 2 > 0. 



The solution to the D 4 quiver consists of five regions and was solved in [22]. Assuming 
that pi > p 2 > p 3 > p 4 > 0, it was found that 

JL_ 1 2jh + 3p, - p 3 1 1 

/2 



1 2pi + 3p 2 - p 3 

4pi (P1+P2) 2 2(pi +p 2 +Ps -Pa) 2(p 4 + p 2 + p 3 + p 4 ) ' 



(3.2) 



We now discuss the solution to the D 5 quiver, consisting of seven regions. We assume 
that k$ > k$ > k 4 > k 3 > k 2 > with k\ = —(k 2 + k 3 + k 4 + 2fcs + 2^) implying 
Pi > P2 > P3 > P4 > P5 > 0. The solution is summarized in Table 1 and Fig. 2 shows 
the eigenvalue distributions and density (further details are given in Appendix B). From the 
information given in Table 1 and (2.3), we find 

1 1 1 

~/2 ~~ 



+ 



2k 2 + 5k 3 + 3k 4 + 4k 5 + 6k 6 2k 2 + 3k 3 + 5k 4 + Ak b + 6/c 6 

4(fc 3 + k 4 + 3k 6 - 2kj) 
(2k 2 + 3k 3 + 3k 4 + 4k 5 + 6k 6 ) 2 

1 1 



9(2k 2 + k 3 + k 4 + 2k 5 + 2k 6 ) 2k 2 + 3k 3 + 3k 4 + Qk 5 + 6k 6 
which, using the relations in (2.7) gives 



18pi 2(pi + 2p 2 ] 



+ 



(2gi + 2p 2 + 3p 3 - p 4 ) 
(Pi +P2+P3) 2 



2(pi + P2 + P3 + Pa ~ Ps) 2(pi + p 2 + p 3 + p 4 + p 5 ) ' 



(3.3) 



Similarly, solving the Dq matrix model as described above leads to a total of nine regions 
and integrating the eigenvalue density gives 



A8 Pl 6(pi + 3p 2 ) 2(p! + p 2 + 2p 3 \ 
1 



+ 



2(pi+p 2 +p 3 ) + 3p 4 -P5 

(Pi +P2 +P3 +P4) 2 
1 



2(pi + + P 3 + p 4 + p 5 - p 6 ) 2(pi + p 2 + P3 + P4 + P5 + Pe) 
for pi > p 2 > ... > P6 > 0. 



(3.4) 
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3.2 General Solution and Polygon Area 



By comparing (3.2), (3.3) and (3.4), we propose that the free energy for D n quivers is deter- 
mined by: 



n—3 n v-^n— 3 . o 

, 2 L&=1 Pb + 3Pn-2 ~ Pn-l 



1 1 

^ = 2 S TS Pb + (n-a- l)p a ' (E^^)' 



+ =J— I , (3.5) 



2 VELlPfc-Pn ELlPfe 

with c„ = 7 ^ and p-i > » 2 > ••• > Pr? > 0. We have verified that this is correct for 

a (n— a— l)(n— a— 2) — ^ A — — 

the D7, D10 matrix models. 

For A-quivers, it was shown in [19] that Vol(Y) can be interpreted as the area of a certain 
polygon. By rewriting (3.5) in a more suggestive form, we will show that there is a certain 
polygon (or rather a cone) whose area is related to Vol(Y) for .D-quivers as well. This 
construction will be particularly useful in Sections 5 and 6. We start by observing that the 
denominators appearing in (3.5) can be written as 



Q-q = ^ (|Pa ~ Pb\ + \Pa + Pb\) ~ 4 |j> a | | a = 1, 71, 
6=1 

n 

a = 2(n - 2) , a n +i = 2}^ \p b \ . (3.6) 

6=1 

The first step in rewriting (3.5) is to combine consecutive terms to get 

Vol( y ) = l( 1 + Pa ~ Pa+l + Pn \ ,g ^ 

Vol(5 7 ) 2 ya ai ^ a a a a+1 o n o n +\)' 

where we have used the relation (2.6). The next step is to introduce the vectors (3 a = (l,p a ) 
together with (3q = (0,1) and (3 n+ i = (1,0). Defining the wedge product (a, b) A (c,d) = 
(ad — be), we can write all the cx a 's in (3.6) in terms of 7 a> & = |/3 a A (3b\ as follows 



<?a = ^ ^ a > b + 7a '- fe ) ~~ 4 7a,n+l ] O = 0, 71 + 1 , (3.8) 
6=1 

where we have also defined /3_ a = (1, — p a ). This finally leads to 

V0l(Y) 1 A 7 a,a+l 



(3.9) 



Vol(,S 7 ) 2^a a <7 a+1 
Now, let us consider the vectors /3 a , a = 0, n + 1 as defining a set of vertices v a given by 



v. 



a = V + ELo A> 
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where Vq is a base point (undetermined for the moment). This set of vertices v a in turn 
defines a new set of edges by the equations v a Ax = 1/2. Then, the set of intersection points 
of consecutive edges, given by w a = (3 a /(2v a A v a +i), together with the origin defines a cone 
C whose area is given by 



1 n 

Area(C) = -^ 



Ah-i A /3 a 



a=0 



v a A v a+ i){v a+ i A v. 



(3.10) 



a+2 



The denominators v a A v a+1 = v a A (v ffl + /3 a ) = u Q A (3 a = (v + J2t=o A /3 a depend on 
the choice of base point vq. Choosing vq = (— n + 2, —1), we can set (v a A v a+ i) = —l/2a a 
leading to 

Vol(Y) 



Area(C) 



(3.11) 



Vol(^ 7 ) 

We also note that by rescaling the cone C by a factor 2/x, we can actually interpret p(x) as 
the height of the cone. In Fig. 3 we show the rescaled cone corresponding to the D§ quiver. 
The x coordinates of the vertices of this cone correspond to the location of the kinks in p(x) 
in Fig. 2. Thus, 1/2 = J dx p(x) = Ap 2 Area(C), from where (3.11) follows immediately. 




Figure 3: Schematic cone for the Ds-quiver. The height of the cone gives the density p(x) in 
the regions defined by the x coordinates of the vertices w a . 



This construction is analogous to the polygon for the A-quiver [19]. The vectors a in 
that case correspond to the (l,q a ) charges of five-branes involved in the brane description 
of the theory. The addition of the two extra vectors /3 and /3 n+ i in the present case seems 
to suggest that one should also include (0,1) and (1,0) branes in the description of these 
theories. 

We would like to comment that solving the matrix model under a different ordering of the 
p's amounts to permuting them correspondingly in the expression (3.9). Moreover, regardless 
of the sign and ordering of p's, the denominators appearing in the expression for Vol(Y) are 
always given by the <r's in (3.8). In the next section we will propose a general expression, 
which is valid for any value of the CS levels and is explicitly invariant under Seiberg duality. 
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4 General Formula for D n Quivers 



It was shown in [23, 24] that the free energy is invariant under a generalized Seiberg duality 
[25, 26]. For ADE quivers, Seiberg duality can be reinterpreted as the action of the Weyl 
group, which acts by permuting and changing the sign of an even number of p's in the case of 
.D-quivers. Thus, we would like to have an expression for Vol(y) that does not assume any 
particular relation among CS levels and is explicitly invariant under Seiberg duality. It was 
proposed in [22] that this can be written as a rational function whose numerator is given by 
J2iz+ det(ai...a n ) 2 YTl=i \ a a • p\, where 1Z + denotes all n-subsets of positive roots. Note that 
under Weyl transformations the a a 's defined in (3.8) are simply shuffled among each other. 
Based on this, we propose that the general expression for the volume corresponding to D n 
quivers is given by 

Vol(y) _ det(ai...a n ) 2 Ua=i K ; p\ , . 

Vol(^)- 2lEo 1 a a ' 

As we will prove below, (4.1) reduces to (3.9) when the CS levels are ordered. 

We recall that in the corresponding formula for A-quivers, the numerator could be inter- 
preted as a sum over tree graphs [15]. In a similar way, we will show now that the numerator 
of (4.1) can be interpreted as a sum over certain graphs known as signed graphs [27] (see also 
[28, 29, 30] and references therein). A graph V = (V,E) consists of a set of vertices V and 
a set E of unordered pairs from V (the edges). A signed graph (r, a) is a graph T with a 
signing a : E — > {+1, — 1} associated to each edge. With these definitions, we can associate 
a signed graph to each term in the numerator of (4.1). Recall that the roots a a for D n are 
of the form (e» ± e,-), where e« are the canonical unit vectors of dimension n and i ^ j. To 
a root of the type (e* — e,) we associate a positive edge (a = 1) connecting the nodes i and 
j in the graph, and to a root of the type (e« + ej) we associate a negative edge (er = — 1). 
Then, we think of the matrix / = {a\...a n ) as an incidence matrix for a diagram with n 
vertices and n edges 3 . Due to Euler's theorem, such graphs must contain loops. If the graph 
contains more than one loop then it must be disconnected. Loops are naturally associated 
a sign as well, given by the product of the signs of all the edges forming the loop. As we 
shall explain below, the determinant in (4.1) selects diagrams containing only negative loops. 
Some examples of diagrams contributing to the numerator for are shown in Fig. 4, where 
dashed lines represent negative edges and solid lines positive ones. 

To understand why the determinant vanishes for diagrams with positive loops, it is useful 
to introduce the operation acting on graphs called 'switching'. Switching is defined with 
respect to a vertex v G V, and it acts by reversing the signs of all the edges connected to 
that vertex. This operation preserves the value of (det I) 2 since it corresponds to multiplying 
some rows and columns of the incidence matrix I by — 1. It is easy to see that by various 
switching operations one can turn any loop with an even number of negative edges into a 

3 Note that due to the absence of roots of the form 2ej, one should not consider edges starting and ending 
on the same node. 
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1 f 



i- 




Figure 4: Some signed graphs contributing to the numerator for D 4 . The first diagram, for 
example, contributes a term 4| (pi + p 2 )(p2 — Pz){Pz — V-±)(Pi ~ Pi) I- 

loop made entirely of positive edges. Then, det / will vanish simply because the columns in 
/ associated to these edges add up to zero. On the other hand, if there are an odd number of 
negative edges in the loop, the above argument does not apply. In fact, one can easily check 
that (det/) 2 = 4 for each negative loop. Thus, we can also write (4.1) as 



V0l(Y) J2( V ,E,a)eT- U(a,b)€E \P* ~ a Pb\ 



Vo\{S 7 



2n 



n+1 - 
a=0 ° a 



(4.2) 



where T denotes the set of signed diagrams with n vertices and n edges (connected or 
disconnected) and no positive loops, L_ is the number of negative loops in the diagram, and 
a the sign of the corresponding edge. Using a generalized matrix-tree formula, we now show 
that (4.2) in fact reduces to (3.9) for p a > p a +i- 

4.1 Generalized Matrix-tree Formula 



We define the n x n adjacency matrix A for a signed graph by: 



n-l 



Ka = ^(lafi + 7a,- 



A a b = -Ja,b + la,-b ■ 



6=1 



The generalized matrix-tree formula [29, 30] states that 

det A = J2 (V!E)(7)eT - 4 L - U( a ,b)eE \Pa - °Pb\ ■ 
By row and column operations we can bring A into the tri-diagonal form: 

/a 1 + a 2 + 2712 —a 2 

-a 2 0-2 + 0-3 + 2723 -03 

-a 3 ; 



A 



— Cn-l 

—a n 



(4.3) 



\ 



1 (~ 



<?n + Pn+l)/ 
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Using the fact that the determinant of tri-diagonal matrices satisfies a recursion relation, we 
have 



detA= - (a n + cr n+ i) det A n _ x - a 2 n det A n _ 2 , 
det A a = (a a + a a +i + 27 a>a+ i) det A a -i - o\ det A a ^ 2 



(4.4) 
(4.5) 



where A a denotes the a x a sub-matrix of A for a = 1, ...,n — 1. Then, using the identities: 
o a +i ~ e a = -2(n - 2 - a)7a,a+i and 



o-l 

E 

d=0 



(n - 2 - rf)7rf |rf+ i _ _J_ 
&dOd+i 2a a 



we can show that the recursion relation (4.5) is solved by 



a-l 



det A a _i = ]~[cr 6 



b=0 d=0 



[a - d)-f d4+ i 



(4.6) 



Using (4.6) in (4.4), we have 



n+l 

det a = 2 n &b 

6=0 



a 



n_ \ v-> (n - ^h^rf+i _ 2o- n ^-4 (n - 1 - d)^ d4+l 

. , / /+ , <t , , , Ft , , 



°n+l / °d °"c2+l 



0"n+l 



d=0 



n+l n— 1 



5lI*E 

n-1 

E 



6=0 d=0 

n+i ri- 1 

6=0 |_d=0 

n+l n 



7d,d+l . / 0"n+l ~ 0"n\ (n - 2 - rf)7d,d+l 



On+1 



<7d<?d+l 



ld,d+i , 1 47„,„+i 1 



c^^d+i 2 cr n+ i 2cr r , 



7d,d+l 



6=0 rf=0^^«H-l 

Finally, substituting (4.3) in (4.2) leads to 

Vol(E) det A 



n 

E 



1d,d+l 



recovering the expression (3.9). 



5 Flavored D n Quivers and the F-theorem 

The F-Theorem [8] states that the free energy (1.1) decreases along RG flows and is stationary 
at the RG fixed points of any three-dimensional field theory (supersymmetric or not). Thus, 
F gives a good measure of the number of degrees of freedom, in analogy with the c-function 
in two dimensions and the anomaly coefficient, a in four dimensions. This theorem was first 
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tested in a variety of field theories [31, 32, 33] and recently it has been proven in [34, 35] for 
any three-dimensional field theory by relating F to the entanglement entropy of a disk-like 
region. Here we check that it holds for the the class of theories we have discussed. We trigger 
the RG flow by adding massive non-chiral fundamental flavors in the UV. By integrating out 
non-chiral flavor fields, there is no effective shift in the CS levels. Thus, we are interested in 
comparing Fik^np) to Fik^Q). The addition of n F ^ in (2.4) introduces no additional 
complications and the matrix model is solved as explained in section 3. We solved the flavored 
D n matrix model for n — 4, 9 leading us to 

Volg) _ 1 ( 701 A 7a,a + l 

Vo\{S 7 ) 2 yaoi&t + n F ) (a a + n F )(a a+1 + n F ) 

By comparing (5.1) with (3.9), it is clear that Fik^np) > F(ki] 0) verifying that 

F uv > F IR , 

in accordance with the F-theorem. 

In terms of the polygon construction discussed in Section 3.2, adding flavor corresponds 
to adding the vector ftp = (0,n F /2) between ft and fti. Then, (5.1) has the same form as 
(3.9) with b = F,l, ...,n in the definition (3.8). 




6 Unfolding D n to A<i n -^ 

Here we provide a check of the formula (3.9), based on the folding/unfolding trick discussed 
in [36], which relates the free energy of various quiver gauge theories when some CS levels 
are identified. It can be used to change the gauge groups from unitary to orthosymplectic 
without changing the quiver or it can be used to change the quiver without changing the 
type of gauge group. Here we will deal with the latter use, as it relates the free energy of 
-D-quivers to that of ^4-quivers. 

When the external CS levels of a D n quiver are identified, it can be unfolded to an Am-^ 
quiver, as shown in Fig. 5. Each internal node in the D quiver is duplicated to give two 
nodes with the same CS level, while the four external nodes combine to give two nodes with 
doubled CS levels. Each node in the A quiver corresponds to a U (2N) gauge group and the 
condition *Yli a n a k a = is automatically satisfied in the unfolded quiver. Using this, it can 
be shown that in the large N limit, Zp, = \J~Z~a and therefore the free energies are related by 
Fp> = \Fa. Here we verify explicitly this proportionality by comparing the formula (3.9) to 
the corresponding formula for A 2n ~5- 

Let us first look at the formula for the D n quiver when external CS levels are identified, 
i.e., k\ = ki = k and k$ = k^ = k'. Due to the relations in (2.7), this is ensured by 
setting pi = p n = 0. Thus, we need the solution to the matrix model with the ordering 
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-0—-0- ' 



Unfolding 



k 




2k' 



Figure 5: Unfolding D n to Am-h- Each node in the A quiver corresponds to a U(2N) gauge 
group. 



Pi > ••• > Pn-i > Pn > Pi > 0. As mentioned at the end of Section 3, this is given by 
permuting the p's in (3.9) accordingly. Then, setting pi = p n = gives 



Vol(Yj 



D 



Vo\(S~< 



-\2l + Y* 7a,a+l + Pn-1 
2 \°2 ^ O-aO-a+1 (^n-l) 



(6.1; 



Now we wish to compare this expression with the corresponding one for A2 n -5 [19], namely 



VoKYa 
Vo\(S 7 



2n-4 



^E 



7a,a+l 



2 f CT a O"a+l 
a=l 



(6.2) 



where a a = YllZ^ \la - Qb\, la,b = \q a ~ Qb\ and ^1=1^1^ = °- Tlie identification of opposite 
CS levels in the A2 n -4 quiver leads to q a = —q2 n -z-a (see Appendix A for details). Then, we 
assume that qi > ... > g n _ 2 > > g n _i > ... > q2n-4 and 



n-2 



2n-4 



n-2 



<jg = [g a - g&l + |g a - = ^(|g a - + |g a + gb|) ; a = l, ...,n-2. 



6=1 ft=n-l 6=1 

Noting that cr a = a a and g a = p a+ i for a = 1, ...,n — 2, we have 



VoU% 
Vol(S 7 



1 / 7a,a+l 

2 V^^^+i 



7n-2,n-l 



n— 3 

E 



7«,a+i 



+ 



C"n-2 Cn-l 
<?n-2 <?1 



2n-5 

E 



7a,a+l 72n-4,2n-3 



^^d^+l cr 2n-4°"2n-3 



n-2 



~2 



(7 



Pn-1 



^2 (o-n-l) 2 

Thus, comparing (6.3) to (6.1) we have 



(6.3) 



Vol(y D ) = -Vol(Y A ). 



(6.4) 



This relation can also be seen clearly in terms of the areas of the corresponding polygons, 
as shown in Fig. 6 (the cone as defined in Section 3.2 has been doubled along the dotted line 
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for visual clarity). The outer polygon corresponds to the A-quiver with opposite CS levels 
identified and the shaded region on the left represents the polygon corresponding to a general 
.D-quiver; when pi = p n = 0, this shaded region expands to fill half of the outer polygon on 
the right. 




Unfolding 




Figure 6: Polygons associated to the -D4 quiver (shaded region) and A% quiver (outer polygon). 
Upon unfolding, Area("P D ) = 1/2 Area^^)- 



Recalling that the nodes of the unfolded A-quiver correspond to U(2N) gauge groups, we 
verify that 

F D _ N 3 / 2 jVol(Y A ) _ 1 
F~ A ~ (2iY)3/2V VoUYo) ~ 2' 



7 Discussion 



In this paper we have studied three-dimensional D n quiver Chern-Simons matter theories by 
using the localization method of Kapustin, Willet and Yaakov in the large N limit. These 
field theories are believed to be dual to M-theory on AdS± x Y, where Y is a tri-Sasaki 
Einstein manifold. We have explicitly solved the corresponding matrix models for various 
values of n, leading us to conjecture a general expression for the free energy and therefore for 
the volume of the corresponding space Y given in (4.1). We have shown that the numerator 
of this expression can be interpreted as a sum over a class of graphs with edges that carry 
a sign, known as signed graphs. Using a generalized matrix-tree formula, we prove that 
for a particular ordering of CS levels, it can also be interpreted as the area of a certain 
polygon, given by (3.9). When external CS levels in the D n quiver are identified, the area 
of this polygon becomes half the area of the polygon corresponding to the A2 n -h quiver, 
in accordance with the unfolding procedure. We have also studied the addition of massive 
flavor fields, showing that when they are integrated out, the area of the corresponding polygon 
always increases (thereby decreasing F), in accordance with the F-theorem. 
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The relevant tri-Sasaki Einstein space for a D n quiver is the base of the hyperkahler 
cone defined by the quotient H 4n_8 ///f/(l) n ~ 1 x SU(2) n ~ 3 . To the best of our knowledge, 
the volumes of these spaces have not been computed. Thus, (4.1) can be considered as an 
AdS/CFT prediction for these volumes. A possible approach to proving the conjectured 
expression for the free energy would be to find the general solution to the matrix model, 
perhaps in terms of the polygon construction presented above, as it has been done for the 
A-quiver in [19]. Some questions which have not been addressed here are whether there is 
a group theory interpretation of the volume formula and whether its denominator can be 
written in a form that is universal for any ADE quiver. 
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A Roots of A m _i and D n 

Here we give some useful information about the roots for A and D Lie algebras. For A m _i 
we choose the following root basis 

a a = e a -e a+1 , a = 1, ...,m - 1; 6 = -e 1 + e m , 

where e a are canonical unit vectors of dimension m. For D n we choose 

a i = e i -e i+ i ) i = 1, ...,n - 1; a n = e n _i + e n , = -(e 1 + e 2 ), 

where are the unit vectors of dimension n. 



Oil ■■■ 9 OL n 




«2n-5 ••• a n-l «1 tt„_l 

Figure 7: Dynkin diagrams for A 2n -5 and D n . 



In Fig. 7 we show the affine Dynkin diagrams for the A and D Lie algebras along with the 
roots associated with every node. At each node, the CS level is given by a ■ q and a ■ p for 
A and D, respectively. The identification of opposite CS levels in the A2 n -5 quiver imposes 
a a ■ q = a 2 „-4-a • q and hence q a = —q2n~3-a f° r a = 1, ...,n — 2. With these conventions, 
unfolding the D-quiver to the A-quiver relates q a = p a+ %. 
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B S 5 



Here we give the detailed solution of the matrix model for the D 5 quiver gauge theory. As 
discussed in Section 3, there are 7 regions defining a generic solution of this model. To keep 
the notation simple, the second index for the four y's corresponding to the external nodes is 
suppressed. 

Region 1: < x < 3{k2+k3+ £ +2k5+2ke) 



(2&4 — ft 3 — k 4 — 2ko) x {2k 2 — k 3 — k 4 — 2ke) x k 3 x 

2/1 ~ 2/6,2 = "j , 2/2 ~ 2/6,2 = ~, , 2/3 - 2/6,2 = ^~ 

4p 4p 2p 

k 4 x (k 3 + k 4 + 2k 6 )x 

2/4 - 2/6,2 = , 2/5,1 - 2/6,2 = 2/5,2 ~ 2/6,2 = 1^ , 2/6,1 ~ 2/6,2 = . 



Region 2: 3 ^ k2+k3+h ^ +2 k B +2k 6 ) 



< X < 



6fc 2 +9fc3+9fc4+12fc 5 +18£!6 



1 (2& 2 — k 3 — k 4 — 2k G ) x k 3 x 

2/1-2/5,2 = --, 2/2-2/6,2 = ^ , 2/3-2/6,2 = "^- 

1 (2/ci + A;3 + + 2A;e) x 

2/4 - 2/6,2 = , 2/5,1 - 2/6,2 = - . , 

2p 2 4p 

1 (2ki -k 3 -k 4 - 2k 6 ) x 

2/5,2 - 2/6,2 = 77 H -. , 2/6,1 - 2/6,2 = . 

2 4p 



Region 3: 6fc2+9fc3+9fc4 / " +1 2fe 5+ i8fc 6 — x — 3( 2 fc 2 +fc 3 +fc4+2A;5+2fc 6 ) 



^ = 3 ! 

1 (fci — fc 2 ) a; (2A?! — 3fc 3 — fc 4 — 2A; 6 ) a; 

2/1 — 2/5,2 = — - , 2/2- 2/6,2 = -1 z , 2/3- 2/6,2 = -1 "j 

2 2p 4p 

{2k\ — k 3 — 3k 4 — 2ke) x 3 k\x 

2/4-2/6,2 = -! - A , 2/5,1 - 2/6,2 = - ;r , 

4p 2 p 

1 (2k x -k 3 -k 4 - 2k 6 ) x 
2/5,2 - 2/6,2 = -g , 2/6,1 - 2/6,2 = -2 — . 
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Rpp-ion 4- ^ < r < ^ 

xvcgiuii <*. 3 ( - 2A . 2+A , 3+A , 4+2A . 5+2fc6 - ) — — 2fc 2 +3(fc3+fc4+2fc 5 +2fe 6 ) 



1 . 1 (2k 3 + h + k 5 + 2k 6 ) x 
Vi- 2/5,2 = ~2' 2/2-2/6,2 = 0, J/3 -2/6,2 = + ^ 

1 (fc 3 + 2A; 4 + A; 5 + 2A; 6 ) x 1 (h - k 2 ) x 
V*-V** = -2 + 2~p ' W,i-W,2 = -2 ^ , 

1 , (h + k 4 + k 5 + 2k % ) x 

2/5,2-2/6,2 = -^, 2/6,1-2/6,2 = -!+ • 

z p 



Rpp-inn "V ^ < t < ^ 

16 2fe2+3(fc3+fc4+2fe5+2fe 6 ) — — 2fc2+3fe3+5fc 4 +4fc 5 +6fe6 



p — fX + xki ; 

1 n (Sk 3 + k 4 + 2k 6 ) x 
Vi- 2/5,2 = ~2> 2/2-2/6,2=0, y 3 -y 6j2 = — 



(fc 3 + 3£; 4 + 2A; 6 ) a; 1 

2/5,1 — 2/6,2 — ^ ' 

(fc 3 + fc 4 + 2fc 6 ) x 



2/4 - 2/6,2 = 1^ , 2/5,1 - 2/6,2 9 



2/5,2 - 2/6,2 = -o > 2/6,1-2/6,2 



2' 2p 



Rpg-ion 6- ^ < r < ^ 

xvcgiuii u. 2 fc2+3fc3+5fe 4 +4A:5+6fc 6 — x — 2fc2+5fc3+3fc4+4fc 5 +6fc6 



p = — + - (6fci - k 3 - 3k A - 2k 6 ) ; 

1 1 (2£; 3 + £; 6 )x 
2/i -2/5,2 = ~2> 2/2-2/6,2 = 0, 2/3-2/6,2= g + ^ , 

11 11 (k 3 + 2k 6 )x 

2/4-2/6,2=2' 2/5,1-2/6,2 = 2' 2/5,2 -2/6,2 = - j, 2/6,1 ~ 2/6,2 = g H ^ 



Rpg-ion 7- ^ < T < ^ 

16 • 2fe>+5fc 3 +3fc 4 +4A;5+6fc6 — — 2fc2+3fc3+3fc4+4fc5+6fc 6 



p = 2p + x (2ki — k 3 — A; 4 — fc 6 ) ; 

1 1 1 

2/i - 2/5,2 = , 2/2-2/6,2 = 0, y 3 -y 62 = -, y^-y 6t2 = -, 

1 11 /c 6 x 

2/5,1-2/6,2=2' 2/5,2 - 2/6,2 = ? 2/6,1-2/6,2=2 + ^-- 

Finally, the last saturation occurs at the end of this region with y ei = y 62 + 1- 
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C Exceptional Quivers 



We have also solved the matrix models for the exceptional quivers Eq, E 7 and Eg. They consist 
of eleven, seventeen and twenty-nine regions respectively. Here we give the corresponding 
free energies for a particular ordering of the CS levels. In Fig. 8, we show our conventions in 
labeling the nodes. 




k\ k 2 k 3 ki k§ k± k 2 k% k 4 k*, fc 6 k 7 k\ ki k 3 k 4 k§ k§ k 7 



Figure 8: Labeling of Chern-Simons levels for E e , E 7 and E 8 . 

The matrix model for Eq gives: 

2 _2{4k 2 + llk 3 + 8k 4 + 4k 5 + 6k 6 + 4k 7 ) 1 

~ {2k 2 + 5k 3 + 4k 4 + 2k 5 + 3k 6 + 2k 7 ) 2 ~ 42(3k 2 + 6k 3 + 4k 4 + 2k 5 + 5k 6 + k 7 ) 

1 1 
~ 77(13A;2 + I2k 3 + 8k 4 + 4k 5 + 3& 6 + 2k 7 ) ~ 3(3k 2 + 6k 3 + 4k 4 + 2k 5 + 5k 6 + 4k 7 ) 

9 9 
~ 6k 2 + Uk 3 + 13k 4 + Qh + 9k 6 + 6k 7 ~ ll(Qk 2 + Uk 3 + 13k 4 + 12k 5 + 9k 6 + Qk 7 ) ' 

under the assumptions that fc 6 > A; 5 > k 4 > k 3 > k 2 > and k 7 > 3k 2 + 6k 3 + 4k 4 + 2k 5 + 4kg. 
% 

The matrix model for E 7 gives: 

2 _8k 2 + 24k 3 + 42k 4 + 4(8fc 5 + 6k G + 3k 7 + 5k 8 ) 
~ {2k 2 + 6k 3 + 10k 4 + 8k 5 + 6k 6 + 3A; 7 + 5k 8 ) 2 

1 1 

2k 2 + 7k 3 + Wk 4 + 8k 5 + 6k 6 + 3k 7 + 5k 8 2k 2 + 6k 3 + 10fc 4 + 9k 5 + 6k 6 + 3k 7 + 5k 8 

1 

_ 180(2A; 2 + 3k 3 + 4k 4 + 3k 5 + 2k 6 + k 7 + 2k 8 ) 

4 

~ 15(4A; 2 + llk 3 + 2(9fc 4 + 8k 5 + 7k 6 + 6k 7 ) + 9k 8 ) 

27 

~ 7(6k 2 + 17k 3 + 28k 4 + 24k 5 + 20k 6 + 9k 7 + 15k 8 ) 

32 

~ 21(8k 2 + 25k 3 + 42k 4 + 32k 5 + 22k 6 + 12k 7 + 27k 8 ) ' 
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under the assumptions that k 7 > k 6 > k 5 > k 4 > k 3 > k 2 > along with Ak 3 + k 4 > 
2k 5 + k 6 ,k 3 + 2k 4 + fc 5 > k 7 , fc 4 + k 5 > k 6 and 3k 8 > 6k 3 + 12fc 4 + 15k 5 + 10k 6 + 5k 7 . 

E 8 

The solution of the matrix model for E 8 gives: 

2 _ 8k 2 + 24k 3 + A8h + 74A; 5 + 92k e + A8k 7 + 64k 8 + 32k 9 
'Jj? ~ (2k 2 + 6k 3 + 12h + I8k 5 + 23k 6 + 12k 7 + 16k 8 + 8k 9 ) 2 

1 

~ 3150(2fc 2 + 3k 3 + Ak 4 + 5k 5 + 6k 6 + 3k 7 + Ak 8 + 2k 9 ) 

1 

~ 2(k 2 + 3k 3 + Qk A + 9k 5 + 12A; 6 + Qk 7 + 8A; 8 + 4fc 9 ) 

1 

~~ 2k 2 + 6k 3 + 13k A + I8k 5 + 23k 6 + 12k 7 + 16k 8 + 8k 9 

27 

~ 7{6k 2 + 18k 3 + 35h + 52fc 5 + 69fc 6 + 38k 7 + 48k 8 + 2Ak 9 ) 

108 

~ 35(12A;2 + 36k 3 + 70k 4 + 10Ak 5 + 138k 6 + 69k 7 + 103k 8 + A8k 9 ) 

36 

~~ 55(12fc 2 + 36k 3 + 70h + 104fc 5 + 138fc 6 + 69k 7 + 103k 8 + 68k 9 ) 

9 

~ 154(6/c 2 + 17(3^ + 4A; 4 + 5k 5 + 6k 6 + 3k 7 + Ak 8 + 2k 9 )) ' 

assuming that k 6 > k 5 > k 4 > k 3 > k 2 > , k 7 > 3k 4 + 6A; 5 + 4/c 6 , 2£; 4 + 4k 5 + 6k e + 9k 7 > 
k 8 , 2k 3 + 4k 4 + Qk 5 + 8k 6 + Ak 7 + Qk 8 > k 9 and 2k 9 > Qk 3 + 12k 4 + I8k 5 + 24fc 6 + lQk 7 + llk 8 . 
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